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(54) Rapidly adaptable channel equalizer 

(57) A method for use in a digital data transfer apparatus, such as a radio telephone, that includes an adaptive 
equalizer, such as an adaptive filter (32), for receiving, correcting, and outputting data symbols, determines 
correction coefficients (34B) for the adaptive equalizer and includes a first step of, for every data symbol output by 
the adaptive equalizer, maintaining a first matrix R that is a NxN channel correlation matrix and a second matrix P 
that is a Nxl cross correlation vector between a received signal and a transmitted symbol (34A). The method 
includes a second step of, for every k* h symbol output (34C, 34D) by the adaptive equalizer, determining a set of 
correction coefficients <C> for the adaptive equalizer based upon the matrices R and P. A third step of the method 
applies the set of determined correction coefficients to the adaptive equalizer so as to optimally correct a 
subsequently received data symbol. The correction coefficients are calculated with an inverse matrix of the 
channel correlation matrix. As a result, the correction coefficients are, in principle, optimal. 
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Rapidly Adaptable Channel Equalizer 



2268377 



This invention relates to data transmission and reception 
equipment , and in particular to methods and apparatus for 
operating an adaptive channel equalizer in such equipment. 

High baudrate data transmission may generally encounter 
intersymbol interference (ISI), where subsequent symbols 
interfere with or overlap one another, making reliable 
reception difficult or even impossible. The ISI might arise 
from amplitude or phase distortion, or in the case of digital 
radio communication, from multipath propagation. 

The mentioned types of distortion are linear and may hence be 
combatted with a channel equalizer. Typically this equalizer 
is implemented as a tapped delay line, finite impulse 
response filter, where the taps or coefficients are 
adjustable, thus making it possible to adjust the filter 
characteristics. 

There are a number of methods for determining the 
coefficients for a channel equalizer. It is known that an 
optimal equalizer can be obtained from a solution of the 
matrix equation: 

RC ~ P, (1) 

where C is a Nxl vector in which the elements are the 
equalizer coefficients, R is a NxN channel correlation 
matrix, and -P is a Kxl cross-correlation vector of samples 
received and symbols sent. 
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In devices typically employing adaptive equalizers, such as 
data modems, it is found that the value N is required to be 
large. As a result, a direct solution of Equation (1) to 
obtain an optimal equalizer requires a large amount of 
calculation capacity. However, in that the required equalizer 
adaptability is generally minimal, relatively simple 
iterative methods, such as a Least Mean Squared (IMS) 
technique, can be employed for solving for suboptimal 
equalizer coefficients. 

When optimal equalizer coefficients are desired, as opposed 
to suboptimal equalizer coefficients, it has been known to 
solve Equation 1 iteratively or recursively by RLS or Kalman 
algorithms. However, these algorithms, in their basic form, 
suffer from a sensitivity to a representation accuracy 
(number of bits) of the numbers used in the calculation. 
Unfortunately, decreasing the sensitivity also rapidly 
increases the requirement for greater calculation capacity. 

In addition, all of the known iterative methods are carried 
out each time a data symbol is received, otherwise the 
equalizer adaptability is adversely affected. 

It is known that a direct solution of Equation 1 may be 
accomplished by inverting the matrix R .and multiplying the 
matrices R~ l and P, or by solving the equation group formed 
by the matrices, of which the latter is a more efficient 
approach. 

According to a first aspect of the invention there is 
provided a method of determining correction coefficients for 
use in an adaptive equalizer means for receiving, equalising, 
and outputting symbols, the method comprising the steps of 
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maintaining a NxN channel correlation matrix and a Nxl 
? crosscorrelation vector matrix P between a received signal 

and transmitted symbols for substantially every symbol output 
* by the adaptive equalizer means determining a set of 

equalizer coefficients (C) for the adaptive equalizer means 
based upon matrices R and P for every k** 1 symbol output by 
the adaptive equalizer means, where k is greater than one; 
and applying the set of determined correction coefficients to 
the adaptive equalizer means to equalize a subsequently 
received data symbol. ■ 

The advantage of the first aspect of the invention is that 
the method exhibits an optimal adaptability, and is also 
insensitive to a number of bits employed for representing 
numbers used in a calculation. The method also provides a 
reduced calculation requirement, without a corresponding 
reduction in the equalizer performance. 

A further advantage of the first aspect of the invention is 
that the adaptive equalizer coefficients are calculated with 
an inverse matrix of a channel correlation matrix. As a 
result, the coefficients are, in principle, optimal. The 
adaptation occurs such that a channel correlation matrix is 
maintained or updated for each received symbol, while the 
correction coefficients are updated only for every k* 
received symbol. In that the method inverts the correlation, 
matrix, and in that the method also calculates the equalizer 
coefficients with a frequency that is less than a frequency 
at which symbols are received, the calculation loading is 
considerably reduced over the known approaches described 
previously. 
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According to a second aspect of the invention there is 
provided an adaptive equalizer means for receiving, 
equalizing, and outputting data symbols, comprising, detector 
means for detecting the data symbols output by the adaptive 
equalizer means and including means for determining 
correction coefficients for the adaptive equalizer means, the 
correction coefficient determining means including means, 
responsive to substantially every data symbol output by the 
adaptive equalizer means, for maintaining a NxN channel 
correlation matrix R and a Nxl crosscorrelation vector matrix 
P between a received signal and transmitted symbols, the 
correction coefficient determining means further being 
responsive to every k th symbol output by said adaptive 
equalizer means for determining a set of correction 
coefficients (C) for the adaptive equalizer means based upon 
the NxN channel correlation matrix R and the Nxl cross 
correlation vector P where k is equal to or greater than two, 
the detector means having an output coupled to an input of 
the adaptive equalizer means for applying the set of 
determined correction coefficients to the adaptive equalizer 
means so as to optimally correct a subsequently received data 
symbol. 

The advantage of the second aspect of the invention is that 
the adaptive equaliser means exhibits an optimal 
adaptability, and is also insensitive to a number of bits 
employed for representing numbers used in a calculation. The 
method also provides a reduced calculation requirement , 
without a corresponding reduction in the equalizer 
performance. 

An embodiment of the invention will now be described by way 
of example only, and with reference to the accompanying 
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drawings in which: 

Fig. 1 is block diagram of a digital radio telephone that 
includes a channel coder/decoder comprising an embodiment of 
the invention; 

Fig. 2 is a block diagram of the decoder portion of the 
channel coder/decoder of Fig. 1; 

Fig. 3 is a more detailed block diagram of the channel 
coder/ decoder, illustrating an adaptive filter and a detector 
used with correlation matrix and coefficient determination 
blocks; and 

Fig. 4a illustrates the N x N channel correlation matrix R, 
while Fig. 4b illustrates a N x N matrix U that maintains the 
matrices S and D. 

There are a number of known methods for solving Equation 1, 
such as Gauss's substitution, Seidel's "and Jordan's iteration 
methods, and Crouch's partition. Another method for solving 
an equation group, attributed to Cholesky and described in 
detail below, is considered to be well suited for use in a 
radio telephone application, and is thus a presently 
preferred technique. It should be realized, however, that the 
teaching of the invention is not limited for use with only 
radio telephone equipment, nor is the teaching of the 
invention limited for use with only the solution approach 
attributed to Cholesky. The Cholesky method is described in 
the text "Digital Communication" by John G. Proakis, second 
edition, published by McGraw-Hill, pages 692-693. 

Referring to Fig. 1, there is shown an overall siaplified 
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block diagram of a digital radio telephone 10. In the ensuing 
description, reference is made only to the reception path of 
the signal, although it should be realized that approximately 
the same signal processing occurs in the transmission path, 
but in a reverse order. 

An antenna 1 provides an input to a radio frequency (RF) 
block 12 which down-converts a received RF signal from an 
assigned frequency, in one or more steps, to a frequency low 
enough {preferably baseband) to be digitally processed. In 
addition, the RF block 12 also accomplishes necessary channel 
filtering. 

A baseband analog signal processing block (RASP) 14 includes 
an automatic gain circuit (AGC) 15, and performs symbol 
synchronization and baseband filtering. The output of the 
HASP 14 is applied to an analog to digital converter (A/D) 
16. 

A channel coder/decoder (CCD) 18 receives the digital 
representation of the synchronized and filtered signal and 
determines channel correction coefficients for an adaptive 
equalizer, as will be described in detail below. 

A speech coder/decoder (SCD) 20 converts the bits received 
from the CCD 18 into a u-law coded Pulse Code Modulated (PCM) 
signal, which is then expanded and converted to an analog 
signal by an audio block (AUDIO) 22. The output of the AUDIO 
block 22 is a final speech signal that is output by a 
loudspeaker 24. A microphone 26 is also coupled to the AUDIO 
block 22 and provides an analog representation of a user's 
speech thereto. The user's speech is then passed through the 
Transmit path, in a manner opposite to that described above, 
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to be eventually transmitted from the antenna 1. 

A control processor (CP) 28 is employed to assign the RF 12 
and the CCD 18 to the selected channel, to manage the 
required protocol between the mobile station* 10 and a base 
"station (not shown), and to manage a required user interface. 
The CP 28 includes or is coupled to a memory 28a that stores 
instructions and data, including the various matrices 
described in detail below. 

Figure 2 is block diagram that illustrates in greater detail 
the constituent components of the CCD 18. It should be 
realized that in some embodiments of the invention more or 
less than the number of components shown in Fig. 2 may be 
employed. 

The input signal from the A/D 16 is lowpass filtered by block 
LPF 30 and is then applied to an adaptive filter (AF) 32. At 
the output of the AF 32 is a detector (DETECT) 34, wherein a 
most likely symbol from the signal constellation and the 
symbol's corresponding bit(s) are detected. As employed 
herein, a symbol is considered to be a signal amplitude and 
phase that is selected from a modulation signal 
constellation. The detected symbol, together with an 
unequalized signal, is employed to update the coefficients of 
the AF 32. The detected bits are deinterleaved (Block 36) and 
convolutionally decoded (Block 38) before being passed to the 
SCD 20 of Figure 1. The operation of Blocks 36 and 38 is 
conventional, and will not be described in further detail. 

The receiver furthermore generally includes a support block 
having an automatic gain control (AGC) , automatic frequency 
control (AFC) • and symbol synchronization circuits that 
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operate in a conventional manner. The output of the symbol 
synchronization circuit is employed to vary the sampling rate 
of the A/D 16. 

Figure 3 is a block diagram showing in greatfer detail the 
Detector 34 and the AF 32. In addition/ there is shown a 
Matrix Calculation block 34A, which operates for every 
symbol, and a Coefficient Calculation block 34B, which is 
switchably coupled to an output of the Matrix Calculation 
block 34A through a switch (SW) 34C. 

It should be noted that all described blocks may be 
implemented with a digital signal processor (DSP), and not 
with discrete circuit elements. For a DSP embodiment the 
Switch 34C is understood to be a logical, as opposed to a 
physical, switching device. 

SW 34C is driven by a counter 34D so as to momentarily close 
once every k symbols so as to couple the output of the Matrix 
Calculation block 34A to the input of the Coefficient 
Calculation block 34B. As a result, the AF 32 coefficients 
are updated once every k symbols, while the Matrix 
calculations are executed for every symbol. The result is a 
nearly optimal AF 32 coefficient update with a -significantly 
reduced computational load, as compared to the conventional 
techniques described above. 

The various calculations performed by the Matrix Calculation 
block 34A and the Coefficient Calculation block 34B are now 
described in detail. 

In the method according to the invention at least once per 
symbol a -sampled signal is received, after which the matrices 
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R and P are maintained by the Matrix Calculation block 34A, 
once per symbol/ according to the following equations: 

R(n+1) - r x R(n) + X(n) x X(n) T {2) 



P(n+1) = r x P(n) + d(n) x X(n)/ (3) 

where X is a vector formed by the samples received* The term 
d(n) is the symbol transmitted from a remote transmitter/ 
typically a base station. The value of the transmitted 
symbol is either known bef orehand, such as during an AF 32 
training phase when a synchronizing symbol preamble is 
received/ or is a detected symbol received during a tracking 
phase. The term r has a value in the range of 0 < r < 1, n 
is a time index # and T indicates a transpose of a matrix. 

It is noted that in a data transfer application R is 
symmetric/ that is R=R H , where H indicates transpose and 
complex conjugate. Thus, the entire R matrix need not be 
maintained/ but only the upper or lower triangle matrix 
portion thereof/ in addition to the diagonal. As a result/ 
the number of elements of the matrix R that must be 
maintained is considerably smaller than the original number 
of elements. 

After updating the matrices R and P according to formulas 2 
and 3/ the following matrix equation is solved by the 
Coefficient Calculation block 34B, 

R x C = P/ (4) 
preferably with Cholesky's method/ described below/ in 
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order to obtain new equalizer coefficients, Cholesky's 
method is presently preferred in that it is a non-iterative, 
non-recursive technique and, thus, reduces the sensitivity to 
the number of bits employed for representing numbers. 

In that the latter operation {Equation 4) requires 
significantly more calculation capacity than the former 
operation (Equations 2 and 3), and because the former 
operation collects all of the existing information of the 
received signal; the latter operation is accomplished less 
frequently than a frequency at which symbols are received. 
That is, the latter operation is performed once every k 
symbols, where k > 2. 

In this manner, the equalizer coefficients that are applied 
to the AF 32 are always optimal when new coefficients are 
calculated for every k th symbol, and the intermediate k, 
k=l.., symbols are corrected with a suboptimal equalizer. 
However, the equalizer adaptability is not adversely 
affected, so long as all the symbols d(n) can be accurately 
identified. 

The value of k is preferably chosen such that symbol 
identification succeeds with a sufficient probability so that 
the audio signal is not noticeably degraded. As was stated, k 
is equal to or greater than two. For some applications it. 
may be preferred to permanently fix or hardwire the value of 
k to some value that is equal to two or greater. 

Alternately, the value of k may be varied during operation by 
the CP 28 over a data path 28A to counter 34D. The value of 
k is made a function of a detected bit error rate (BER) so. as 
to maintain k at a value that provides an acceptable BER for 
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a given set of transmission and reception conditions. 

Cholesky's method for solving a group of equations, which is 
a presently preferred technique for implementing the 
Coefficient Calculation block 34B, is as follows. 

STEP A: 

As seen in Fig. 4a, the matrix R is partitioned into three 
matrices S, D, and S H so that 

R «= S x D x S H , (5) 

where S is a triangular matrix of R, such as the lower 
triangular matrix, and D is a diagonal matrix of R. The 
partition is carried out according to the following formulas: 

D{1) = R(l,l) (6) 

S(i,j) = (R(i,j) - 2 S(i,k)D(k)S(j,k))/D(j), (7) 

k-1 

where 1 < j < i-1 and 2 < i < N, and 

D(i)=R(M)-2s(i,l)S(i,i|D(l!) ) (8) 

M. 

where 2 < i < N 

STEP B: — • • ... 

next, the matrix Y-D-x-tf 1 x C is defined and the following 
group of equations: 
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s x y = p. 

is solved with the following formulas: 



(9) 



Y(l) = P(l), and ' (10) 

i-l 

Y(i) = P(i) - 2 S(i,j)Y(j), (11) 
j-l 

where 2 < i < N. 
STEP C: 



Next, the set, of equalizer coefficients C are determined from 
the following group of equations: 

S H x C - D~ l x Y (12) 

with the following formulas: 

C(N) = Y (N) / D(N) , and (13) 

C(i) = Y(i)/D(i) - XS(j,i)C(j), (14) 

j-i+i 

where 1 < i < N-l. 

For minimizing a memory requirement, and for facilitating 
memory addressing, a NxN matrix U may be defined, as shown in 
Fig. 4b, in which the matrices S and D are maintained in the 
following form: 

U = S x D7D + D" 1 + S 8 -!, (15) 
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where I is a diagonal unit matrix. This matrix storage 
'1 scheme enables in-place factorization. It is noted that the 

matrix Y can physically occupy the same memory locations as 
j the matrix C so as to further minimize memory 28a 

requirements . 

While the invention has been particularly shown and described 
with respect to a preferred embodiment thereof, it will be 
understood by those skilled in the art that changes in form 
and details may be made therein without departing from the 
scope and spirit of „the invention. 

The scope of the present disclosure includes any novel 
feature or combination of features disclosed therein either 
explicitly or implicitly or any generalisation thereof 
irrespective of whether or not it relates to the claimed 
invention or mitigates any or all of the problems addressed 
by the present invention. The applicant hereby gives notice 
that new claims may be formulated to such features during 
prosecution of this application or of any such further 
application derived therefrom. 
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Claims 

1. A method of determining correction coefficients for use 
in an adaptive equalizer means for receiving, equalising, and 
outputting symbols, the method comprising thfe steps of: 

maintaining a NxN channel correlation matrix and a Nxl 
crosscorrelation vector matrix P between a received signal 
and transmitted symbols for substantially every symbol output 
by the adaptive equalizer means; 

determining a set of equalizer coefficients (C) for the 
adaptive equalizer means based upon matrices R and P for 
every k th symbol output by the adaptive equalizer means, 
where k is greater than one; and 

applying the set of determined correction coefficients 
to the adaptive equalizer means to equalize a subsequently 
received data symbol. 

2. A method according to claim 1, wherein the step of 
maintaining includes the steps of operating on the matrix R 
in accordance with the expressions: 

R(n+1) = r x R(n) + X(n) x X(n) T , and 

P(n+1) - r x P(n) + d{n) x X(n), 

where X is a vector formed by received symbol samples, d(n) 
is a symbol, the value of which is either known beforehand or 
is detected, r has a value in the range of 0 < r < 1, n is a 
time index, and T indicates a transpose of a matrix. 
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3. A method according to claim 1 or claim 2, wherein only a 
v diagonal matrix D and either an upper triangular matrix or a 

lower triangular matrix of the matrix R is maintained. 

4 . A method according to any preceding claim, wherein the 
. step of determining includes a step of: 

solving the expression 

RC » P 

by the steps of, . 
STEP A: 

partitioning the matrix R into three matrices S, D, and 
S H so that 

R = S x D x S H , 

where S is a triangular matrix of R and D is the 
diagonal matrix of R, the step of partitioning including a 
step of solving the formulas: 

D(l) m R(l,l) 

M 

S(i,j) = (R(i,j) - 2S{i,k)D(k)S(j,k))/D(j), 

k-l 

where 1 < j < i-1 and 2 < i < N, and 
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M 

D(i) = R(i,i) - 2S(i,k)S(i,k)D(k), ^ 

k-l 

where 2 < i < N; 
STEP B: 

defining a matrix Y - D x S H x C, and solving the group 
of equations: 

S x Y = P 

in accordance with the formulas: 

Y{1) = P{1), and 

Y(i) = P(i) - IS(i,j)Y(j), 
j-i 

where 2 < i < N; and 
STEP C: 

determining the set of correction coefficients C in 
accordance with: 

S H x C = D" 1 x Y 

by the use of the formulas: 

C(N) « Y (N) / D(N), and 

- . : . . k 

C(i) - Y(i)/D(i) - S S(j,i)C{j), 

j-in 

-where Ki< H-l, 
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5, A method according to claim 4, and further including a 
step of defining a NxN matrix U for maintaining the matrices 
S and D as: 

U = S x D-D + D'^-I, 
where I is a diagonal unit matrix. 

6, An adaptive equalizer means for receiving, equalizing, 
and outputting data symbols, comprising: 

detector means for detecting the data symbols output by 
the adaptive equalizer means and including means for 
determining correction coefficients for the adaptive 
equalizer means; 

the correction coefficient determining means including 
means, responsive to substantially every data symbol output 
by the adaptive equalizer means, for maintaining a NxN 
channel correlation matrix R and a Nxl crosscorrelation 
vector matrix P between a received signal and transmitted 
symbols, 

the correction coefficient determining means further 
being responsive to every k** symbol output by said adaptive, 
equalizer means for determining a set of correction 
coefficients (C) for the adaptive equalizer means based upon 
the NxN channel correlation matrix R and the Nxl cross 
correlation vector -P where k is equal to or greater than two, 

the detector means having an output coupled to an input 
of the adaptive equalizer means for applying the set of 



18 



determined correction coefficients to the adaptive equalizer 
means so as to optimally correct a subsequently received data 
symbol . 

7. An adaptive equaliser means according tt> claim 6, 
wherein the correction coefficient determination means is 
operable on the matrix R in accordance with the expressions: 

R(n+1) a r x R(n) + X(n) x X(n) T , and 

*(n+l} = r x P(n) + d(n) x X(n), 

where X is a vector formed by received data symbol samples, 
d(n) is a data symbol sent, the value of which is either 
toown beforehand oris detected, r has a value in the range 
of 0 < r < l, n is a time index, and T indicates a transpose 
of a matrix. 

8. An adaptive equalizer means according to claim 6 or 
claim 7, and further including memory means for storing a 
diagonal matrix D and either an upper triangular matrix or a 
lower triangular matrix of the matrix R.. 

9. An adaptive equalizer means according to any of claims 6 - 
to 8, wherein the correction coefficient determination means 
solves an expression RC = P by a first operation of: 

partitioning the matrix R into three matrices S, D, and 
S" so that 

* - S x V x S H , 
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where S is a triangular matrix of R and D is the 
diagonal matrix of R, the operation of partitioning including 
the operations of solving the formulas: 

D{1) = R(l,l) 

. 5-1 

S(i,j) = (R(i,j) - 2 Sli^JDIkjSlj.kJJ/Dlj), 

k-l 

where 1 < j < i-1 and 2 < i < N, and 

j-i 

D(i) = R(i,"i) - 2 S^kJS^JclDfk), 
where 2 < i < N; 

by an operation of defining a matrix Y = D x S H x C, and 
solving the group of equations: 

S x Y ■ p 
in accordance with the formulas: 

Y(l) = P(l), and 
i-i 

Y(i) = P(i) - 2 S(i,j)Y(j), 
j-i 

where 2 < i < N; and 

by an operation of determining the set of correction 
-coefficients C in accordance with: 

£ H x C « D" 1 x Y 
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by the use of the formulas: 
C(N). = Y(N) / D(N), and 

C(i) = Y(i)/D(i) - 2S(j,i)C(jJ, 

j-i+l 

where 1 < i < n-1. 

10. An adaptive equalizer means according to claim 9, and 
further including memory means for storing a NxN matrix U so 
as to maintain the matrices S and D as: 

U = S x D-D + D"*+S*-I, 
where I is a diagonal unit matrix. 

11. An adaptive equalizer means according to any of claims 6 
to 10 comprising an adaptive filter. 

12. A radio telephone comprising an adaptive equalizer means 
operable an accordance with a method according to claims 1 to 
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according to any of claims 6 to 10 



A radio telephone comprising an adaptive equalizer means 



14. A radio telephone according to claim 12, further 
comprising: 

means for receiving and for transmitting an RP sional 
that is moduXated with digital data represenLSvTof Sech; 
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and 

coding and decoding means for applying and extracting 
j the digital data as data symbols/ said coding and decoding 

means including, for a receive path. 

15. A method for operating an adaptive equalizer means 
substantially as hereinbefore described and with reference to 
the drawings. 

16. An adaptive equalizer means substantially as 
hereinbefore described and with "reference to the drawings. 
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